Fluctuation Effects in Steric Reaction-Diffusion Systems 
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We propose a simple model for reaction-diffusion systems with orientational constraints on the 
reactivity of particles, and map it onto a field theory with upper critical dimension d c — 2. To 
two-loop level the long-time particle density N(t) is given by the same universal expression as 
for a nonsteric system, with N(t) ~ t~ d ^ 2 for d < 2. For slow rotations of the particles we find an 
intermediate regime where N(t) ~ t~ d ' A , with the crossover to the nonsteric asymptotics determined 
by the rates of rotations and reactions. Consequences for experiments are discussed. 
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Diffusion-controlled reactions are ubiquitous in Nature and span a wide range of phenomena — from chemical 
catalysis in cell biology to matter-antimatter annihilation in the early Universe. In these processes the time to for 
a particle to diffuse to its nearest neighbor is much larger than the characteristic time tn for the particles to react, 
once in proximity. As a result, large density fluctuations may develop at late times, resulting in anomalous kinetics 
with universal scaling behavior jl| . 

Our understanding of diffusion-controlled reactions has advanced significantly in the last few years. Results obtained 
by new theoretical approaches (like renormalization group methods or mappings to integrable models H ) can now 
be tested against high-quality experimental data, including measurements on low-dimensional processes where effects 
from fluctuations are most pronounced However, little attention has been paid to the modeling of fluctuation 
effects in steric reaction-diffusion systems, i.e. systems where the reactivity between particles depend on their relative 
orientation. Examples include diffusion-controlled enzyme reactions |^|, effective reaction-diffusion models of protein 
folding H , as well as generic solution kinetics of molecules with nonuniform surface reactivity Q . Here the chemically 
active part of a reactant is located in a specific region, usually represented by an active spot on the sphere modeling 
the reactant. Two particles react only if their active spots touch, signaling a successful "docking" for reaction. Despite 
the importance of this class of processes, the theoretical effort so far has been limited to estimating how the reaction 
rates effectively decrease due to orientational constraints [Qj8| . The lack of a more comprehensive theory reflects the 
subtlety of the problem, which here makes a faithful treatment quite difficult. 

In this Letter we take a new route and propose a minimal model for steric reaction-diffusion systems, specifically 
annihilation processes. This model - in spirit an "Ising-type" model for steric reaction-diffusion - has the advantage 
that it can be analyzed by field-theoretic renormalization group (RG) techniques, allowing for precise predictions 
about the long-time kinetics. In spite of its simplicity the model contains some very interesting information. In 
particular, for sufficiently slow internal rotations of the particles our results reveal that steric reactivity in dimensions 
d < 2 drives a crossover from a density decay N(t) ~ t~ d / 4 at large intermediate times to N(t) ~ t~ d / 2 at later times 
(with logarithmic corrections in d = 2). 

To define the model we partition configuration space into cells Ar, and time into intervals At. The particles are 
represented by spheres with an active spot that is either "up" or "down" with respect to a fixed reference direction, 
and "A" ("B") denotes a sphere with the spot up (down) g|. Thus, the orientation of a sphere is restricted to two 
states, A and B. The number of A and B particles in the cells at a time step nAt, n = 1,2, are given by the 
sets of cell occupation numbers {a}„At and {6} n At respectively. We assume that the microscopic dynamics lead to 
a Markov process po| with waiting time interval At, and assign probabilities that within At a particle in a given 
cell will react (if the cell is nonempty), jump to a neighboring cell with or without rotation (with rates D tra nsi+rot 
and D respectively), or stay in the cell with rotation (with rate D rot ) or without. In the continuum limit D rot and 
Dtmns+rot fold into one effective parameter which we refer to as D rot , with the continuum limit translational diffusion 
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constant denoted by D. At first disregarding the influence of microscopic rotations we assume that two particles in a 
cell may react and annihilate with some rate 5q only if one is A and the other B, with 5o an average over all possible 
active ( "spot-to-spot" ) and nonactive encounters within the time interval At. Since microscopic rotations in the time 
interval At may turn a nonactive encounter into an active one we also assign a nonzero reaction rate Ao to the cases 
where both particles initially are of A or B type. The assignment of values to the rate constants is immaterial, but 
clearly Ao will approach 5o for sufficiently fast microscopic rotations. Thus, to "optimize" effects from steric reactivity 
we assume that Ao <C So, and also D ro t/D <C 1. This limit of "slow" rotations is indeed realized in many experimental 
systems [||. 

The model thus constructed is equivalent to a generalized two-species annihilation system where A and B particles 
diffuse by translations and "flips" A <-» B and react according to jll]] : 

A + A^(l) A + B^® B + B^®. (1) 

Without rotations (Ao = and D rot — 0) one immediately obtains the well-known A + B — > universality class for 
the asymptotic mean density: n a (t) = rib(t) ~ t~ d / 4 jl^fl. To explore the more realistic case including rotations we 
consider the master equation 

j t P(c, t)=J2 Rc^ c P(c', t) - £ R c ^P(c, t). (2) 

c' c' 

Here P(c, t) denotes the probability of a given configuration c = ({a}, {&}) at time t, and R c '^ c is the transition rate 
from state d = ({a}', {&}') into c, determined by D, D rot , S and A . We take P(c, 0) as a Poisson distribution, with 
averages denoted by n a fl and nb t o for A and B particles respectively. In the absence of anisotropics or external fields 
we expect that n a ,o = f^bfi, but for now leave the initial densities unspecified. 

Following Refs. the master equation (2) can be mapped onto a field theory in the continuum limit, with A 

and B particles described by scalar fields a and b respectively. Introducing <fi = (a+b) /2, tp = (a— b)/2 (with conjugated 
fields defined without the factor of 1/2), and corresponding initial densities n^.o = (n a fi+nb,o)/2,n^fi = (n a ,o — "f>,o)/2, 
we obtain the action 

S = J d d x J dt[4>{d t - V 2 )0 + $(dt - V 2 )V> + + (A + (5 O )00 2 + (Ao - S )W 2 
+2AoV^ + + + + V 

+\ a 4>Tp(j)ij}\- J d d x(n<f,fi(i){x, 0) +7ty,oV>(a:, 0)), (3) 

where we have rescaled the time variable, Dt — > t, and where p ~ D rot /D. Using (|^), the mean densities can be 
calculated as n c (t) = (c(t))s (c = a, b), where ( )s denotes an average with respect to the action S. The rate 
constants, here playing the role of coupling constants, become dimensionless when d = 2, suggesting this as the upper 
critical dimension above which mean field theory ("normal reaction kinetics" jjj) should be valid, a result confirmed 
by standard power counting techniques [jfE]]. Hence, we concentrate on d < 2. These dimensions can be explored by 
using RG improved perturbation theory in d c = 2 and then continuing the results to d < 2 via an e-expansion. As 
seen in Eq. (^), rotational diffusion acts as a mass perturbation in the "0-sector, and for simplicity we here focus on 
the "strong" steric limit where this term is omitted, commenting on the general case as we go along. 

Using dimensional regularization, divergences in perturbation theory show up as poles in e = 2 — d. Analogous 
to the single-species case A + A — > |^6|, these poles can be removed by simply renormalizing the rate constants: 
Ao — > Afl and 5q — > Sr. The renormalized rate constants are defined by primitively divergent vertex functions 
evaluated at an arbitrary momentum scale k, which is further used to define dimensionless bare and renormalized 
rates go = K~ t 5o, ho = K~ e Xo and gn = K~ e 5n, ha = k~ il \r, respectively. These are related by the exact result 

9o , h Q 

l +9o/g l + ho/h* 

with g* = h* — r(e/2) _1 (87r) d / 2 the fixed point to which both rate constants flow. We are interested in the large-time 
mean densities n a (t) and rib(t) which depend in principle on the initial densities n 0) o and Ubfl, and the bare rates Ao 
and So, expressed through hn,gn and k. However, the densities cannot depend on the arbitrarily chosen scale k, and 
this implies a Callan-Symanzik equation 

n c (t; g R , h R ; n afi , n bfi ; k) = (n 2 ty d/2 n c (K~ 2 ; gR{t),h R {t);n afi {t),n bfl {t); k) , c = a,b (5) 



2 



with running coupling constants 



~ 9R(t) - i + mr^ hR{t) - irmF^ 2 ' (6) 

and densities n Cj o(i) = (n 2 t) d ^ 2 n c _ , with c = a, b. The crossover times are given by t* s = (g*/So) 2 ^ e and t* x = (/i*/Ao) 2//e 
respectively. The fact that the rate constants flow to the same fixed point coupling g* = h* has some unexpected 
consequences, and dramatically simplifies calculations as compared to the two-species case A + B — > with Ao = 
— a notoriously difficult problem within the RG [p[ . 

To calculate the mean densities n a (t) and n&(i), the perturbation series has to be organized in such a way as to 
make it well-behaved when the initial running densities flow to infinity in Eq. (|5|). At tree level this can be done by 
simply "dressing" the initial densities n c o (c — a, b) 

x = x+ S + ... (7) 

and free propagators Gq(*i — W,P) = ®{t\ — t%)e~ p2 tt 1 ^ 2 } 



In these diagrams a thin (thick) line corresponds to a free (dressed) propagator, with a cross representing an initial 
density. The dressed tree-level densities (a(t)) and (b(t)) (with averages taken only over tree-level diagrams) are 
solutions of the coupled classical rate equations 

|(a(t)) = -(Xo(a(t)) 2 + S {a(t)){b(t))), ±(b(t)) = -(X (b(t)) 2 + S {b(t))(a(t))), (9) 
which for equal initial densities n a = n^ g = n yield the prediction 

^ = ^ = n n\x^ - ( 10 ) 

1 + n (A + do)t 

This result is drastically changed by fluctuations, as revealed by a loop expansion of the densities. To carry it out 
we need the dressed tree-level propagators, which for equal initial densities can be calculated exactly. In the <p-tp 
basis: 

G,( P M - t 2 ) = 6(U - t 2 )e-^-^ ( ±± n f° ± S ;\ t2 ) Ai , £ = ^ (ID 

\1 + n (X + 5 Q )tiJ 

with A<p = 2 and = 2Ao/(Ao + So). At one-loop level there are two diagrams which contribute to (4>) (with ( ) 
now denoting a full average): 

x x 
(0)i- ioop = -^^^ + " (12) 

with a full (dashed) line denoting a <j)(tp) propagator. After performing the momentum loop integrations and taking 
no — ► oo, one obtains 

Mi-loop = \(^ty d/2 I dxx 1 -^ 2 [ aly[y 2 + z 2 y 2z ](l - V y d l 2 (13) 
1 Jo Jo 

with t\ = tx and t% = txy, and with z = (Xr — Sr)/(Xr + 5r). Eqs. ^ and (||) imply that z — » for large times, 
thus killing the second term in the square bracket in Eq. (|l^). By an e-expansion, and using that (ip) = for equal 
initial densities, we thus obtain for large times: 



n a (t) = n b (t) = - 



1 In 8?r - | 
1 

Aire 8tt 



-d/2 



t^oo, A^O, al<2 (14) 
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This result, which we have checked to two-loop level, yields the same universal expression for the total density 
N(t) = n a (t) + njj{t) as for a nonsteric system fl6f . The inclusion of rotational diffusion with p ^ in Eq. (g) 
does not influence the expression in Eq. ([14]): since p does not renormalize, the long-time asymptotics of the theory 
stays the same [0. The fact that the asymptotic density is blind to steric effects and exhibits not only the same 
decay exponent —d/2 as a nonsteric system but also the same amplitude is a striking manifestation of universality in 
reaction-diffusion processes. 

Eq. ( |f3| ) reveals how Ao controls the universality class of the model. In addition to the pole in 1/e, there is a 
second potential divergence coming from the y-integration due to the flow no — > oo under the RG. However, as no 
flows, z flows as well unless X = 0, and the asymptotic value of z determines the fate of the integral: for z > —1/2 
(z < —1/2) the y— integral is finite (divergent). If Ao — one has the A + B — > universality class 0. For this case 
z = — 1 < — 1/2 for all times and the integral diverges. On the other hand, for any A ^ 0, z — ► and accordingly 
z > —1/2 for large times, producing a well-behaved expansion. This gives the A + A — > universality class of Eq. 
(|l4|) . Our results imply that even if Ao is arbitrarily small, the system will eventually behave as if Ao = 5q ( non-steric 
limit). How long does it take for this to happen? Since z = —1/2 is the critical value which controls the no dependence 
of perturbation theory, and thus the crossover to the A + A — > universality class, it is natural to define a crossover 
time t z such that z(t z ) = —1/2. Eq. fafy together with the definition of z shows that 



Since tg -C t z ~ t* x for small Ao, there is a large intermediate time interval with a strong no- dependence where gn is 
close to its fixed point value g* , while tiR is still small. A strong no-dependence with gn ~ g* implies scaling with the 
exponent — d/A when hn = (no rotations) Qj. Since a small hji cannot change this behavior, it will be present also 
for slow rotational modes with Ao ^ 0. The inclusion of rotational diffusion with p ^= in Eq. (|^) will decrease t z , but 
for sufficiently small p there will still be an intermediate regime with a strong no-dependence, supporting scaling with 
the exponent — d/A jl7fl . We expect this property to hold generically, also for more realistic systems with continuous 
degrees of freedom. By analogy, we also expect a similar crossover from n(t) ~ t~ d l 4 to the usual mean-field behavior 
n(i) ~ t^ 1 for steric systems in dimensions d > 2. Work on this problem is in progress. 

What causes the crossover from the t~ d l 4 decay to the t~ d / 2 deacy for small Ao? The most probable scenario is 
that A and B particles begin to segregate into A and B rich regions, a process during which the A + B — > reaction 
is dominant (—d/A decay exponent). Once A and B regions are formed, however, the A + B — > reaction becomes 
unimportant and the A + A — > and B + B — > reactions become dominant despite the small value of Ao , resulting 
in the —d/2 decay exponent. 

Figure 1 shows the results of a numerical simulation of the master equation @ in d = 1 [jisj . The uppermost curve 
corresponds to A = and S = 10 (A + B — > model with decay exponent —d/A), while the lowest curve corresponds 
to the case with equal reaction rates A = S = 10, representing the asymptotic limit of equal rate constants as given in 
(^) . As predicted, for small values of A (intermediate curves) the decay exponent (slope of the curves) changes from 
— d/A for early times to —d/2 at late times, with crossovers consistent with ( |l5| ) (scaling with 1/A 2 for small A). 

The present model can be generalized to an arbitrary number of orientations N by introducing the "particles" 

A±, .... An, and reactions A a + Ap —4 0. This generalized system can be shown to belong to the same universality 
class as the present N — 2 model as long as all rates \ a ,p are nonzero |L7j. However, if one or several of the rates 
A a ,/3 vanish, a different scaling behavior results (most probably the t~ d / 4 asymptotics of the "A+B" model, although 
this is still to be established). 

The recent breakthrough in preparing and monitoring low-dimensional chemical reactions — using carbon nanotubes 
as effectively one-dimensional "test tubes" ]l9| , p0| — opens up the possibility to experimentally test our predictions. 
The essential elements to be present in a laboratory set-up are that (i) the reactions are sufficiently fast to ensure the 
diffusion-controlled condition t^ << tp, (ii) the reactants have a highly anisotropic surface reactivity with sufficient 
slow rotational diffusion to make the system strongly steric, and (iii) the reactive products precipitate out of the 
system on a short time scale compared to tjj and hence do not influence the kinetics An observation of the 

crossover behavior predicted in the present work would be extremely interesting as it would exhibit, for the first time, 
a reaction-diffusion process supporting two distinct time scales with different anomalous scaling. 
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FIG. 1. Density deacy obtained from a numerical simulation of the master equation in (2) in d = 1 with 10 5 cells. Here 
8 = 10 and n° A — n% = 3.5 particles/site. 
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